Determination of the Carrier-Envelope Phase of Few-Cycle Laser Pulses with 

Terahertz-Emission Spectroscopy 
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The availability of few-cycle optical pulses 
opens a window to physical phenomena occur- 
ring on the attosecond time scale. In order 
to take full advantage of such pulses, it is cru- 
cial to measureii2j2j4 and stabilise^"^ their carrier- 
envelope (CE) phase, i.e., the phase difference be- 
tween the carrier wave and the envelope function. 
We introduce a novel approach to determine the 
CE phase by down-conversion of the laser light to 
the terahertz (THz) frequency range via plasma 
generation in ambient air, an isotropic medium 
where optical rectification (down-conversion) in 
the forward direction is only possible if the inver- 
sion symmetry is broken by electrical or optical 
means®'^-^'*'^--'^*'. We show that few-cycle pulses 
directly produce a spatial charge asymmetry in 
the plasma. The asymmetry, associated with THz 
emission, depends on the CE phase, which allows 
for a determination of the phase by measurement 
of the amplitude and polarity of the THz pulse. 

The ability to measurei*2iM and stabihseii^ the CE 
phase, which is also often named absolute phase and de- 
fined by Eq.^ (Methods), of few-cycle laser pulses, i.e., 
the ability to completely control ultrashort high-intensity 
electromagnetic fields, provides access to a whole series 
of novel physical phenomena, and creates options for 
their coherent control. An example is the manipula- 
tion of the electron recombination process during "re- 
collision"—, with the aim to fully control the generation 
of both intense VUV higher-harmonic radiationiSii^ as 
well as of high-intensity single attosecond pulses. The 
progress has led to the establishment of a new field 
of research, "attosecond science" , with highlights of re- 
cent work including the efficient tomography of molecu- 
lar wave functions by coherent diffraction of rescattered 
electrons--" as well as the preparation of high-intensity, 
monochromatic femtosecond electron beams evolving 
from "bubble generation"—. Many of the phenomena 
studied depend critically on the CE phase. Exam- 
ples include the dynamics of above-threshold ionisation 
(ATI)'', light-induced non-sequential double ionisationi^, 
Gouy phase observationii, and quantum interference in 



photocurrentsi^. An exciting prospect is the application 
of the CE phase as a coherent-control tool, e.g., to steer 
bound electrons in moleculesi^. 

The state-of-the-art method for CE phase determi- 
nation of amplified laser pulses is stereo ATlSLSi. In 
this letter, we present a new approach based on down- 
conversion into the THz frequency range in a laser- 
generated ambient air-plasma. The concept relies on the 
detection of the THz electromagnetic pulses emitted dur- 
ing photo-ionization of the air by the focussed few-cycle 
pulses. The method is easy to apply, because it does not 
require a vacuum chamber and relies on standard opto- 
electronic THz detection technologyi^. 

THz-pulse emission from laser-generated plasmas has 
been investigated in the past, but only with much longer 
laser pulses. There, the transient electric currents re- 
sponsible for the THz radiation have their origin either 
in ponderomotive forcesS^, the Coulomb force of an ex- 
ternally applied DC bias^, or an optical "AC bias" 
achieved by superimposing the light pulses with their own 
second-harmonic radiatiori£i2i2ii£. The THz- wave genera- 
tion mechanism for the latter has been discussed in terms 
of a four-wave-mixing (FWM) processi^iiS. For the 
present case where the pulse spectra are octave-spanning 
FWM results in a low-frequency (fl 0) polarisation: 

P {Q = uji + UJ2 - uJs) oc 1 (uji) I (w2) I (^s) COS {(Pce), 
because the CE phases (pcE of all three input waves are 
the same (as all frequency components arise from the 
same near-transform-limited few-cycle pulse). The CE 
phase is hence directly imparted onto the single-cycle 
THz pulse and can thus be determined via electroop- 
tic measurement of the THz field amplitude. The same 
argument concerning the transfer of the CE phase to the 
THz field can be readily generalised to the case of higher- 
order non- linear contributions. 

Whilst such a wave-mixing process illustrates the mech- 
anism for CE phase-sensitive THz generation, we will 
give an alternative description based on static-tunnelling 
theory^"^, which explains the origin of non- linear polarisa- 
tion on a microscopic level. For simplicity, we assume the 
air to consist only of N2. Figure ^ displays the calculated 
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FIG. 1: lonisation rate for vs. electric field strength based 
on the static-tunnelling theory as a function of the optical 
field strength, respectively the light intensity. The calculation 
was performed with the help of Eq.0of Sec.|I] An lonisation 



rate of 10' 



corresponds to a probability of 100 % for a 



molecule to be ionised during a time interval of 1 fs. Note, 
that the probability for lonisation of a molecule has a value 
of 1% per fs at a threshold field of S.STO'" V/m (vertical 
black line), and that the lonisation probability per fs already 
becomes 100% at S.S-IO^" V/m, i.e., at less than twice the 
threshold field. 
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lonisation rate Wtun for N2 at a pressure of 1 atm, respec- 
tively the lonisation probability rate of a single molecule. 
The strongly non-linear intensity dependence of the lon- 
isation rate translates into temporal lonisation patterns 
as illustrated in Fig. El The upper panels display the two 
pulse shapes considered here, both with the same pulse 
duration of r = 6fs (FWHM), the same centre wave- 
length of A = 800 nm, and the same peak electric field of 
the envelope of — 4.2T0^"V/m, but with different 
CE phases, (pcE = and (pcE = 7r/2. The middle panels 
of Fig.|21exhibit the time evolution of the lonisation rates 
for the two cases. A common signature is that only the 
three or four most intense half-cycles of the pulses ionise 
a significant amount of molecules. The ionisation rate at 
each half-cycle, however, depends sensitively on the CE 
phase. This becomes even clearer when we calculate the 
temporal evolution of the plasma density. We distinguish 
two quantities, p'^{t) and p~{t), which describe the den- 
sities of free electrons produced at positive and negative 
polarity of the optical laser field, respectively. The total 
density of free electrons is given by ptotai = (t) + (t) ■ 
The lower pair of panels of Fig.[5]shows that for (fcE — 0, 
the total plasma density generated at negative field po- 
larity is significantly smaller than the density for positive 
polarity. For tpcE = 7''/2, the densities only differ little 
after the pulse has passed through the focus. 

Electrons generated at positive field polarity experi- 
ence a displacement in the positive field direction due 
to acceleration in the laser field and incoherent scat- 
tering at neighbouring molecules within a few fs. Cor- 



FIG. 2: Optical electric field, ionisation rate and plasma den- 
sities as a function of time for CE phases ifcE ~ and 
ipcE = 7r/2. The plasma densities and p~ , generated at 
positive respectively negative field polarities, are calculated 
with Eqs.|5|and|21of Sec.|I] The ionisation threshold is the ar- 
bitrarily chosen electric-field value where the ionisation prob- 
ability is 1% per fs. 



respondingly, electrons generated at negative polarity 
are displaced in the negative direction. Such a collec- 
tive charge displacement after the passage of the pulse 
(limt_>oo (p^ (0 — p {t)) 7^ 0), creates a quasi-DC macro- 
scopic dipolc moment. Because of its transient character, 
it acts as the source of an electromagnetic pulse, whose 
polarity and field amplitude are determined by the ratio 
i?^ = limt^ 00 {{p'^ — P~)/{p'^ + p~)), aquantity which is 
proportional to the maximal generated dipole moment. 
A value of i?^ = ±l corresponds to the highest possi- 
ble THz pulse amplitude, with the polarity of the signal 
depending on the sign of i?^. For =0, THz emis- 
sion vanishes. Figure|31 displays i?^ as a function of tpcE 
for laser pulses with Gaussian envelope of various dura- 
tions T. One finds that the magnitude of R'^ is directly 
related to the CE phase with a nearly cosine-like depen- 
dence. The contrast, however, rapidly vanishes when the 
pulse duration rises towards 10 fs. 

Whilst the theory for stereo ATI requires modelling of 
both the quiver motion and recoUision of electrons, this 
is not necessary for THz emission from a plasma gen- 
erated at atmospheric pressure. The reason for this is 
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FIG. 3: Asymmetry parameter of the plasma density vs. 
CE phase for several pulse durations. is proportional to 
the macroscopic quasi-DC dipole moment which is responsible 
for the emission of the THz radiation. The pulse intensity 
is chosen such that a plasma density ptotai of 10% of the 
molecule density Pgasis obtained. Inset: Maximal value of 
at LpcE = as a function of the pulse duration. 
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that the coherent motion of the majority of electrons is 
terminated by ultrafast scattering before recoUision and 
the THz signal is proportional to the asymmetry in the 
tunnelling process. Therefore the signal in our method 
is conceptually and in contrast to stereo ATI directly re- 
lated to the CE phase in a way such that a cosine-like 
pulse {}pcE = 0) will produce a maximum asymmetry in 
ionisation rates and also a maximum THz-signal. The 
introduced THz method would measure directly the ab- 
solute CE phase, if the plasma would be a point object. 
However, the plasma is an extended object— and we 
expect a CE phase variation by dispersion during pulse 
propagation and a concomitant interference of the THz 
signals generated along the beam path. The strength of 
this effect for the moment is unknown, but will be inves- 
tigated in future measurements and simulations. 

An experimental test of these theoretical findings was 
performed by THz-emission spectroscopy with phase- 
stabilised 8-fs amplified laser pulses (see in Methods a 
detailed description of the set-up displayed in Fig.0J. 
We measured the THz field value from a laser-generated 
air plasma at a fixed temporal position (that with the 
largest amplitude of the THz transients) while continu- 
ously varying the CE phase by slowly ramping the set- 
point of the control loop (with a slope of 2n per minute). 
An exemplary time series spanning 8 minutes is shown 
in Fig.[S^. The CE phase data measured in the fre- 
quency domain by two f-to-2f interferometers^ (CEP2 
and CEP3) are compared to the THz data. The THz sig- 
nal varies with the ramped CE phase, as expected from 
the model calculations shown in Fig.|3| 

In a second experiment, the CE phase was kept con- 
stant in the oscillator and laser amplifier (CEP2), but 
varied with the help of a piezo-driven glass-wedge pair 
behind the pulse compressor and in front of both the 



FIG. 4: Sketch of the Experimental set-up. a) Schematic of 
the laser system. AOM, acousto-optical modulator; CEPl, f- 
to-2f interferometer used for CE phase stabilisation with fast 
feedback; CPA, chirped pulse amplification and compression; 
CEP2, f-to-2f interferometer used for compensation of slow 
drifts of the CE phase; gas-filled fibre used for spectral broad- 
ening; CEP3, f-to-2f interferometer only used for CE phase 
measurements; b) THz system employed for electro-optical 
THz detection. Inset: A temporal THz wave-form obtained 
with applied DC bias. 



third f-to-2f interferometer (CEP3) and the THz system 
(also with a phase ramping of 27r per minute). The re- 
sults are depicted in Fig.[Sl3. The measurements again 
consolidate the expected change of the THz signal with 
varying CE phase. 

In order to provide a more quantitative analysis of the 
measurements, we show in Fig.|3; the THz signal as a 
function of the CEP3 data (open dots), averaged over a 
74-minute time series, where we ramped the phase of the 
oscillator as in Fig.[S^. The error bar indicates the stan- 
dard deviation of the THz-signal values. The correlation 
between THz signal and CEP3 compares well with the 
dependence predicted in Fig.|3| for 8-fs pulses (full line 
in Fig.O;). The slight deviation of the theoretical data 
from the experimental ones may be a consequence of a 
non-Gaussian shape of the envelope. 

The status of the THz approach for CE phase determi- 
nation can be considered as a proof-of-principle. Never- 
theless, to classify its current precision, we compare it to 
stereo ATI, which achieves a 3tT error of lOOmrad for 10 
seconds of data acquisition and for 6.5-fs 20- /i J pulses—. 
The precision reached with the THz system is TOOmrad 
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FIG. 5: a) Variation of tlie CE phase of 8-fs pulses obtained 
by slowly ramping the set-point of the control loop. The 
CE phase is measured by a f-to-2f interferometer (CEP2) 
directly behind the laser amplifier, another f-to-2f interfer- 
ometer (CEP3) behind the recompressor, and by the THz 
measurement system via electro-optical detection of the THz 
radiation generated in an ambient air plasma. The integra- 
tion time constant of the lock-in amplifier in the THz mea- 
surements was 100 ms. In addition, a sliding average was 
calculated numerically, bringing the effective integration time 
constant to 5 sec. Dispersion effects result in an arbitrary 
offset between the CE phases measured by CEP2, CEP3 and 
the THz system, b) Variation of the CE phase of 8-fs pulses 
by a change of the dispersion in the beam path, realised by 
two variable glass wedges, c) Adjacent average of the correla- 
tion of the CEP3 and the THz data. The window chosen for 
averaging was ±0.5 radians of CEP3. For the averaging, the 
total data pool recorded within 74 min. is taken into account. 
The error bar depicts the standard deviation of this averag- 
ing. The solid line shows the results of model calculations of 
R'^ for an 8-fs pulse (constant scaling factor between theory 
and experiment). The peak electric field used for calculations 



is E%t = 4 • 
the experiment. 



10 V/m corresponding to the expected value in 



for 10 seconds of data acquisition for 8-fs 66-/xJ pulses. 
Considering the factor-of-4 improvement to be achieved 
with 6.5-fs instead of 8-fs pulses (see Fig.OIl, the THz ap- 
proach compares well with the much more mature ATI 
method already now. The sensitivity of low-repetition- 
rate THz-detection can be enhanced significantly by box- 
car techniques or by THz interferometry with intensity 
detection over a large THz-photon bandwidth. 

In summary, we report (i) frequency down-conversion 
in an isotropic medium without breaking the symmetry 
of the non-linear medium by external means, and (ii) CE 
phase measurement of sub-lO-fs amplified optical pulses, 
both by detection of the THz radiation emitted from a 
laser-generated ambient-air plasma. The theory of tun- 
nelling ionisation explains the findings well, in particu- 
lar the down-conversion (THz generation) by a spatially 
asymmetric process of ionisation and charge acceleration. 
THz-emission spectroscopy possesses the potential to be- 



come a commonly employed tool for determination and 
long term stabilisation of the CE phase for amplified few- 
cycle pulses. 



I. METHODS 

Ionisation model employed in the simulations 

We base the quantitative description of the photo- 
ionisation process of atoms and molecules on static- 
tunnelling theory-3., with the following expression for the 
ionisation rate: 



Wt'u 




X exp 



(1) 



where uja — 4.13-10^^ 1/s is the atomic frequency unit, 
Ea = 5.145-10^^ V/m the atomic unit of the electric field, 
and Eopt the optical field strength. J7j^„ = 13.6 eV and 
^ion = 15.6 eV denote the ionisation potentials of atomic 
hydrogen and molecular nitrogen, respectively. The sim- 
ple tunnelling model is known to describe the experimen- 
tal data for N2 well if the intensity of the laser pulses is 
above 10^"* W/cm^, and to achieve a fair agreement down 
to 5-1013 W/cm"^^. 

Temporal evolution of the ionisation 
In order to determine the temporal evolution of the 
plasma current, we compute the plasma densities p~^(t) 
and p~{t) 
equations: 



by solving the following coupled differential 



dt 



dt 



P'{t) 



Wtun{Eopt{t)) ■ {plas " P+ (*) " 

yiEopt{t) > 



HEoptit) <0, 



Wtun{Eopt{t)) ■ {plas 



(2) 



P+{t)-p-{t)), 



0, 



HEoptit) <0 
'^iEopt{t) > 0. 



Here, 



rgas 



(3) 



2.7- 10^^ cm~'^ is the density of N2 molecules 
at atmospheric pressure assuming air to only consist of 
^2- Eopt{t) denotes the electric field strength of the laser 
pulse. We neglect plasma recombination which occurs 
on a much longer time scala^^. We assume the following 
temporal wave-form of the laser pulse: 



Eopt (t) 



^opt 



Ait) cos I —t + ificE 



(4) 



where A{t) is a Gaussian envelope. Whilst such a product 
separation of carrier and envelope can become ill-defined 
as the pulse duration approaches the sub-cycle range^ 
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it is perfectly reasonable for the shortest pulse durations 
we consider here. 

Light source 
The experiments were performed with a light source as 
sketched in Fig.0^. A commercial available Ti:Saphire 
Multipass CPA System "Femtopower Compact Pro" 
(Femtolasers, Vienna) is seeded with a self-build laser 
oscillator (pulse energy 4 nJ, pulse duration 20fs, centre 
wavelength 810 nm, CE phase-stabilised with f-to-2f self 
referencingi (CEPl)) and driven at a repetition rate of 
2.8 kHz. Its output (500 ^J, 30fs, @805 nm) is focused 
into a gas-filled hollow fibre (inner diameter 200 /im, 
length 82 cm, 4.5 bar Neon). A small fraction of the CPA 
output (few /iJ) are used for a f-to-2f interferometeu^ 
(CEP2) to correct for slow CE phase drifts in the am- 
plifier. The fibre output is then compressed to its band- 
width limited duration of 8 fs by the use of dispersive 
mirrors and a prism-compressor sequence with two fused 
silica prisms, with final output parameters (83 /xJ, 8 fs, 
@750 nm (centre wavelength)). 

THz-emission spectroscopy system 
80% of the final output is fed into the THz detection set- 



up, which is sketched in Fig.^J). 10% (the residual 10% 
are unused) is used for a third f-to-2f interferometer— 
(CEP3), which is tracing CE phase changes for direct 
comparison with the THz data in respect to an con- 
stant arbitrary offset. The THz system is a typical THz 
pump-probe measurement system with balanced electro- 
optical THz detection^'^, the whole set-up is fitting on a 
30x45-cm^ breadboard. The pump beam is mechanically 
chopped at half the laser repetition rate. The plasma is 
generated in the focus of an off-axis paraboloidal mirror 
with an effective focal length of 76 mm. For optimi- 
sation of the alignment and a first determination of the 
zero-delay point for the THz detection, an additional DC 
bias field can be applied to the focal region^. The probe 
beam passes along a variable delay line and is focused 
collinearly with the THz beam into the electro-optical 
ZnTe detector crystal. The THz wave-form as measured 
with a DC-biased air plasma is depicted in the lower left 
part of Fig.^lD. For the CE phase measurements, the 
delay stage of the probe beam is positioned at the zero- 
delay point where the strongest THz signal is obtained, 
and the external bias field is switched off. 
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